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In this paper we suggest a consistent approach to derivation of generalized Fokker-Planck equation (GFPE) 
for Gaussian non-Markovian processes with stationary increments. This approach allows us to construct the 
probability density function (PDF) without a need to solve the GFPE. We employ our method to obtain the 

lO ■ GFPE and PDFs for free generalized Brownian motion and the one in harmonic potential for the case of 

power-law correlation function of the noise. We prove the fact that the considered systems may be described 
with Einstein-Smoluchowski equation at high viscosity levels and long times. We also compare the results with 

>-Q I those obtained by other authors. At last, we calculate PDF of thermodynamical work in the stochastic system 

^ I which consists of a particle embedded in a harmonic potential moving with constant velocity, and check the 

work fluctuation theorem for such a system. 
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1. Introduction 



i-rt ■ The theory of Markovian Brownian motion is successfully used in describing a great variety of 

P3 ' experiments and observations [l|-y|. However, it remains an inapplicable model for the majority 

O , of natural systems where a characteristic time of thermal fluctuations is comparable to that of a 

Brownian particle (Gaussian non-Markovian systems), or where the processes are strongly non- 
Gaussian (either Markovian or non-Markovian), all of which results in the fact that the long-time 
mean squared displacement does not grow linearly in time any more, (^x^{t)j ex i^. This phe- 
nomenon is called anomalous diffusion, namely, when /i < 1, the system is said to be subdiffusive, 
and when /i > 1, it is superdiffusive. Evidently, when /i = 1 we have an ordinary Brownian motion. 

^^ , There are two paradigmatic models describing anomalous diffusion: continuous time random 

l^ ' walk (CTRW) and fractional Brownian motion (FBM). The CTRW approach was developed by 

Montroll and Weiss in 1965 [5| for a description of the electric charge transport in a disordered 

f^ ' medium (amorphous semiconductor) [6|. This model considers the independent identically dis- 

tributed couples of random space-time steps whose PDFs belong to the domain of attraction of 
Levy stable laws. 

Recently, the Markovian Levyprocesses in external fields were studied by means of Langevin 

L^ ■ and fractional kinetics technique [7|-|9|. 

''^ , The second model (FBM) was introduced by Kolmogorov in 1940 [10| and later studied by 

C^ ' Yaglom [11]. The name "fractional Brownian motion" belongs to Mandelbrot and van Ness who 

suggested a stochastic integral representation of this process [U] . FBM is a continuous centered 
non-Markovian Gaussian process X^^\t) with covariance function 

(x(^)(t)X(^)(f')) = D {f" + t'^" -\t~ t'|2«) , (1.1) 

or, at large times, 

'x'^"\tf) = 2Dt'", (1.2) 
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where H is Hurst index, < iJ < 1, and D is the generaUzed diffusion coefficient of the dimen- 
sion [D] — cm^/sec^^. Previously, the problem of particle escape from the potential well in the 
framework of this model was considered in paper [13| by using the method of numerical simulation 
of Langevin equation with fractional Gaussian noise Y''^\t). The latter is a non-Markovian sta- 
tionary random process which is defined as the time derivative of FBM and whose autocorrelation 
function exhibits a slow decay at infinity as (F^^^ (i)r(^) (0)) « 2DH{2H- l)t'^"-^, in contrast to 
white Gaussian noise, where (Y^-^/^^t)Y^^/^\t')'^ = 2DS{t — t'). Power spectral density for white 
noise does not depend on frequency, otherwise the noise is called a coloured noise. 

The pioneer work of deriving a differential equation (in essence a Fokker-Planck equation, FPE) 
describing ordinary Brownian motion (OEM) was done by Lord Rayleigh [lj|, within the approach 
of an absence of external potential and an overdamped discrete motion of a heavy Brownian 
particle. A more consistent method was developed by Fokker, Smoluchowski and Planck (a detailed 
historical sketch may be found in '4!|). However, when dealing with the coloured noise case, the 
above-mentioned approaches are no longer valid. 

The most common example of derivation of one-dimensional Fokker-Planck equation for colou- 
red noise may be found in paper [4]; for one-dimensional case it was done in [15]; for a particular 
case of a linear oscillator it was obtained and studied in [16i] . The multi-dimensional case was 
considered in [17|. 

The theory of generalized Brownian motion (GBM) finds its applications in many problems 
of modern physics, biophysics and astronomy. Indeed, polymers [18-20], elastic chains and mem- 



branes [19|, I2ll424| and rough surfaces |25l427j can be described by a continuum elastic model which 
accounts for their general stochastic behavior; it was recently shown that the probe particle in such 
systems performs FBM [28|, |29| . The fluctuations of magnetic field in the turbulent plasma of the 
Earth's magnetospheric tail turn out to have colour: in the range of frequencies lu ^ 10^^ Hz they 
have the properties of flicker-noise (their power spectrum is proportional to l/w). When uj is about 
10~^ Hz, they are a brown noise with the tendency of "blackening" at lower frequencies, see, e.g., 
the paper [30] and works cited therein. Moreover, a similar situation is known from experiments 
in laboratory plasmas: it was found that the power spectra of the saturation current, electrostatic 
potential fluctuations, and the turbulence-induced flux measured in various plasma devices [31] 
have power-law dependencies. At high frequencies, an asymptotic power fall-off of the fluctuation 
spectra with characteristic decay indices close to 2 was denoted; at intermediate frequencies, the 
decay indices were about 1, gaining a weak frequency dependence at the lowest frequencies. 

Another important application comes from single-molecule dynamics. In paper [32j it is shown, 
that the experimental data of the distance fluctuations between the two components of fluorescein- 
tyrosine complex can be described within the framework of the Langevin equation with harmonic 
potential and coloured source possessing correlation function (CF), which decays as t~o'5i±oo'i'_ 

Below we present a consistent method of derivation of a multi-dimensional generalized Fokker- 
Planck equation for linear stochastic systems driven by coloured Gaussian noise paying special 
attention to the case of coloured Gaussian noise with power-law correlation function. 

2. Basics of the method 

We use the approach to obtaining an ordinary Fokker-Planck equation for linear systems with 
delta-correlated noise described in monograph |33,] as the basis of the suggested method for deriva- 
tion of the generalized Fokker-Planck equation. The paper continues and extends the previous 
studies [3J] where we considered the GFPE for exponential and power-law correlation function 
restricting ourselves only to space-homogeneous case. Here we study a more general problem for 
the power-law correlation function. For the integrity and clarity of presentation, we give a full 
description of the method, as well. 

First, let us write Langevin equations in multi-dimensional form: 

i^ = -a^k^k+Y,{t)+K,. (2.1) 

Here S,i is the generalized coordinate, aik is the coefficient matrix, Yi is the external noise, Ki is 
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the regular constant force; the dot above ^i stands for time derivative. Let the initial conditions 
be ^i (t = 0) = £,i (0). Then, the formal solution of (PTTj) is 

t 

e. (i; C (0)) = (e--^*)^^. C, (0) + fdr (e^f*-^)) ^ . {Y, (r) + K,) , (2.2) 



where a = \\aik\\ is matrix composed from the elements of aik . The probability density function 
(PDF) of the value £^i in the moment of time t with the fixed f (0) is evidently a multi-dimensional 
Dirac delta-function: 

/ (?, t; e (0)) = <5 (e - C (i, e (0))) ^Y[SiC,-^, (t, C (0))) . (2.3) 

i 

In case we do not know the exact ^ (0), but their initial PDF / (^ (0) , 0), the PDF at an arbitrary 
moment of time t will be of the shape: 

/ it = / de (0) ,/ (C (0) , 0) (<5 (e - e (t, e (0)))) , (2.4) 

where (. . .) stands for J drpy (t) . . . and py (t) is the PDF of noise. By using the n-dimensional 
delta- function representation 6 (^) = (27r) " J dgexp (ig^) and taking into account (|2.2p . we have: 

(^ (e - C (t, e (0)))) - (27r)-" J dqG (g, t) exp [ig (^ - 6""*^ (0))] , (2.5) 

where 

G (q, t)^l exp J ~iq f dre--^^*-") {Y (t) + K)\) . (2.6) 



Here we should remark that we use a matrix notation and the hat indicates that the value is a 
Fourier image. 

Expanding the PDF into Fourier integral 



/(C,i)-(2^)""ydge''«/(9,*), (2-7) 

we get from (EH) and (^3)) : 

f{q,t)^G{q,t)f(^{e--'fq,0), (2.8) 



^'^ I 1C1 O TV! O+T-TV +-T'0 TTOT-VJ^Od/-] ^- /~\ /U '"'■'' 



where (e""*) is a matrix transposed to e" 

Hereinafter we assume the random process Yi (t) to be a stationary Gaussian process, so that 
the following relations are true: 

(^ii (^l) • • -^isn+l (^2n+l)) = 0, 

{Yij^ (ti) . . .Yi^^{t2n)) = 2^ 9iii2 i^l ^ ^2) . ■ . gi2„_ii2„ {t2n-l — hn) , (2.9) 

where the summation is executed by all possible pair compositions of ii , ii ; 12 , ^2 ; ■ ■ • ; *2n , i2n ■ 
The number of such pairs is (2?! — 1)!! = 2n!/n!2". gi-^i^ {ti — ^2) is a certain function of time 
difference. 

Using the exponential function series expansion for (j2.6p and keeping in mind (j2.9p we have: 



6(M) ^ f '-')'•' (^»)' 



dii / di2g.(e-"*0.,'?™K"*0™/3j'(^i-^2) 







(2n)! n!2'' 
xexp j-iq,ydr(e-''^)^^.i^,j . (2.10) 
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Here and below for simplicity we write gji (ti ~ t-2) instead of gj-^i.^ (t\ ~ t-2). 
Introducing the value 

t t 

97J.™ (i) ^\JdhJdt2 {c-^*-)^^ (e~'^*^)„;5,7 ih - h) , (2.11) 







we get for G {q,t): 



G {q, t) = cxp 



-qiQr.i'OJlim (t) - iqi / dr (e ''^)^^ Kj 





(2.12) 



2.1. Fokker-Planck equation 

It may be easily proven that the value G {q, t) obeys the following relation: 



dG dG 

-17- + qiaik-^ h iKiqiG = -qiq„iDirn {t) G {q, t) , 

Ot dqu 



(2.13) 



where 



'^irn {J' ) 



dOT, 



dt 



+ aik'y^km (t) + Umk^ik (t) ■ 



On the other hand, due to an obvious equality 



(2.14) 



d_ 
di' 



^J ( (e--^*)^ q, O) = -qa^f ((e-^Yl^ « 



(2.15) 



and (|2.8p we can conclude that the function / {q, t) obeys the same equation as (|2.13p which after 
the inverse Fourier transform yields: 



J^i TTZ = -77- [aimkmj [C,, t)\ + XJim (t) 



dt 



d^i d^i 



^^^^^r, 



actually being the generalized Fokker-Planck equation (GFPE). 
Now, let us simplify the expressions (|2.14l) . Noticing that 



(2.16) 



a^k[e )^^ 
dgji (ti - ^2) 



i^-n. 



d_ 

dt ^ "y 
dgji (ti - ^2) 
dt. 



and 



dDJli. 



t t 

If... _„,_^ . _„,^ ,. .. 1 



- J dt, (e--0., (e-")™z 5.' ih^t) + -J dt, (e--)^^. (e--^)^, 5,. (t ~ h) . 




we arrive at the expression 



^.m (0 = 2/ dti (e-''*0„-5,m dill) + (e-"*0™,5., (I^il) 



(2.17) 
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2.2. Probability density function 

The advantage of the described method is that there is no need in solving the Fokker-Planck 
equation to obtain the probabihty density function since we have constructed it imphcitly at the 
stage of the GFPE derivation. Indeed, according to (|2.8p , knowing the Fourier image of the initial 
PDF / {q. 0) = S' {/ (^, 0)} we can easily get the PDF for an arbitrary moment of time: 



= i?-i{G(g,t)/((e-'^*)^q,0)} 



(2.18) 



where G {q, t) is given with (I2.12p . 

However, the expressions (I2.1ip may be rather complicated for direct calculations regarding, 
e.g., a power-law correlation function g (ii — ^2)- By means of integration variables change we get 
a much more usable expression, because now the internal integral does not contain the correlation 
function: 



t t-T 





(e-^) 



ml 



K'^"). 



-a{T-\-T) 



ml 



3. Generalized Brownian motion 

Let us now apply the derived formulae to the specific stochastic system: the generalization of 
the classical Brownian motion with the external random force is a stationary Gaussian noise with 
long memory effects. 

3.1. Free generalized Brownian motion. Spatially homogenous case 

First, we investigate a simple system described with the following Langevin equations: 

da; 

d^ = "' 

~lv + Yit), (3.1) 



dv 
di 



where x(t) is particles coordinate, 7 is friction constant, Y{t) is Gaussian external noise with 
{Y{t)Y{t')) = g{t — t'). The coefficient matrix ai,„ is 

-1 

7 

and gij (ii — ^2) = Si2Sj2g (^1 — ^2) • The solution of the homogenous system p.ip yields: 



x{t) = Xo ^ 
v{t) = Woe" 
Comparing these expressions with p.2p we get 

1 




-t7 



-Wo, 



■t7 



1 - e"*^ 



-t7 



(3.2) 



(3.3) 



In what follows we restrict ourselves to the power-law noise correlation function of the form 



5(r) 



|r|^r(l-/3) 



(3.4) 
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with < /3 < 1, which is actually the asymptotics of the CF for fractional Gaussian noise. Note, 
that at /3 — > 1, we get the delta- function limit ^(t) — > cS{t) [35|. 

Now, we can write out the exact values for the coefficients 1)ij and 9Jt,j , see equations (|2.17p 
and (|2.19p . respectively: 



SiiW = 0, 

Si2(i) = S)2i(i) 



ct 



1-/3 



2,r(2-/.)[^ + (^"^)^^(^^)]-^^^'"^^ 



2)22(i) 



^ r(i-/?) ' 



-C7^ 



5Hi2(i) = 

a«22(t) = 



73r(2-/3) 7'r(3-/3) 2^ 

c (-e-^*T + 2e'*T) M(l - ^, 2 - /3, ^7)^^-^ cE^ft7)ti-^ 
^ 273r(2 - /3) ^ 273r(l - (3) 

'.{l-c-*^)t^-^ c (-e-^*^ + e-*^) ti-^M(l - /3, 2 - /3, ^7) 



272r(2 - (3) 2-f^T{2 - /3) 

1 /3_2 cii-^E;3(i7) ce-2*Tii-'3M(l-/3,2-;9,t7) 
" 27r(l - /?) 27r(2 - /?) 



2^7^ 



Here E^(i) is an integral exponential function 



(3.5) 
(3.6) 

(3.7) 



(3.8) 

(3.9) 

(3.10) 



E^t) = / dp 



-tp 



P 



fi 



and M(a, 6,i) is Kummer's confluent hypergeometric function: 

1 



M(a,&,i) = 



T{b) 



T{a)T{b~a) 



due*"u"'i(l-u) 



h-a-l 



(see, e.g. [3a|). ^^ 

The generalized Fokker-Planck equation (|2.16p for this case yields: 

^^ = ^[7-/(«,i)]+C22(0^^;^ 



When 7^ ^ 1, the latter expression takes the form 

df{v,t) d 



dt 



*'rt-P\ d^f{v,t) 



= ^ [7./(M)] + - ^7^ - ^.(^3^ j ^ 



(3.11) 



(3.12) 



According to the procedure described in section \n^ the PDF f{v,t) with the initial condition 

f{v, 0) — nS{v — vq) reads 



fiv,t) 



2^ACT 



exp 



{v-e *^uo) 



where 



(7 = — C7' 

2 ' 



L.,/3-2 _ ct^~'^^pitl) _ ce-2*T^i-/^M(l-/3,2-/3,t7) 
" 27r(l - /3) 27r(2 - /3) 



(3.13) 



(3.14) 
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3.2. Free generalized Brownian motion. Spatially inhomogenous case 

Now we examine the same system but with inhomogenous initial condition /(a;, v, 0) = nS{x — 
xo)5{v — vq). AU the values for S)ij(i) and dJlij{t) clearly, remain the same as in the previous 
section, but the GFPE and the PDF do change: 



df df d[vf) ^^ d^f ^ d^f 

oxov ov^ 



dt dx dv 



Again we construct the solution with the procedure explained in section 



fix,V,t\LU = 0) = 



where 



and 



^^exp ~^^nit)(p^^^^^~voe 



p^ — (l - e *^) - X + xq 

7 



a = OTn(t)OT22W-»ti2(i)' 



2 

v\ - 



4OTii(i) 



(3.15) 

(3.16) 

(3.17) 
(3.18) 



3.3. Generalized Brownian motion of linear oscillator 

Here we study the most general system, though restricting ourselves to the case of a harmonic 
potential U {x) = uj'^x'^ /2 . The pair of Langevin equations now have the following form: 

dx 

d^ ^ "' 

dv 

dt 



= -711 - oj'^x + Y (t) , 



(3.19) 



where x{t) is particles coordinate, 7 is friction constant, u is frequency of the linear oscillator, Y[t) 
is the external noise. The coefficient matrix introduced in equation (j2.ip is 







1 



uj^ 7 
Again, gij (ti — ^2) = <^J2<5j2.g (ti — ^2) • The solution of the homogenous system (j3.19|) yields 

v{t) = Aie-^*/2e"*/2_^^2e-^*/2e-^"/2, 



xit) = -^[v{t)+^vit)] 



2lj2 



A^^-lt/2^nt/2 (^ + O) + A2e-^*/2g-i«/2 (^ _ Q) 



(3.20) 



where Ai and A2 are constants depending on the initial conditions and here we introduce the value 
U = +-^72 — 4^2. Assigning x (0) = xq and v (0) = wq , we get 



Ai = - 



2uP-Xq + Vq (7 - ?L) 



A. 



2a;^xo + f p (7 + f^) 

2rj 



(3.21) 



Now, substituting the latter expressions into equations p.20p and comparing the result with 
the formal solution (|2.2p without the integral term (since we are looking for the solution of the 
homogenous system), we find 



e-at ^ g-7t/2 



, ,'f7A 7 , /m\ 2 , /m 

cosh — H smh — — smh — 

2yr2 \i j n \2 

2uP- fnt\ fnt\ 7 fnt 

— sinh -j cosh - --sinh - 



(3.22) 
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The final step before proceeding to tlie GFPE and the PDF evaluation is to obtain the exact 
expressions for dJlij (t) and the generalized diffusion coefficients J),y (t) {i,j — 1, 2) for our power- 
law correlation function 



A straightforward integration of equation (|2.19p with (|3.4p gives 

^g-2Qpt |-^^^^ _ 2e"^Lj2) M(l - (3,2- /3, apt)t^'^ 



Wtii(i) 



2-fUj^n^T{2 - 13) 

ce-*'' (ape*"7 - 2lu^) M(1 -/3,2-P, a,nt)t^-P 
27w2f72r(2 - /3) 



2-fUj^nT{l - (3) 27w2f]r(l - /3) 2apa„7w2f^ 

^12 W = -^^Hiw:^-^^'''' (e2^-*-e*^)M(l-/3,2-/3,api) 



Wl22(t) 



where 17 = ^72 - 4^2, Cp = (7 + r2)/2, a™ = (7 - fi)/2. 

According to equation (|2.16p . the GFPE for such a system reads 



'-^^ ^ si-"/<-')i-^"'-^- 



-vf{x, V, t)] + — [{u'^x + jv)f{x, V, t)] 



where 



(3.23) 



2f72r(2-/3) 

, /g2api _ gt(7+20)\ j^(-j^ - /?, 2 - ;3, a,„i) , (3.24) 

apce-2'^p* (7 - 2a™e*") M(l - ^, 2 - /3, apt)^!-'' 

27r22r(2 - (3) 

a^ce*("-"°^) (e*"7 ~ 2ap) M(l - /3, 2 ^ /3, a,^^)^^-^ 
27f72r(2 - 13) 

apcEi3{apt)t^-'^ a^cEi:i{a^t)t^-^ (o^ - gg,) c 
27J7r(l-/3) 27f7r(l-/3) 27^} ' ^ ^ 



92 f 32 f 92 J- 

Dn^ + (Dr2+D2r)^ + D22^, (3.26) 



Sii = 0, (3.27) 

^ _ ^ _ a„apc[E^(apt)-E^(a„<)]ii-'3 c[a^(7 + r!)-a(^(7-r!)] 

®^^ - ^''- 2u^nni-f3) + A^JHi ' ^^-^^^ 

^ c [g^Eff (a„t) - flpEff (apt)] t^-^ c (g^ - a^) 

®^^ - f7r(i-/?) + n ■ (^-^^^ 

At this stage we may compare these diffusion coefficients to that obtained in the paper by 
Wang and Masoliver [l^- We consider only the case of the external driving noise (see section 3.2 
of the mentioned paper). To establish a connection with our GFPE and equation (W29) (here 
the letter "W" indicates the reference to the equation from the paper [IB]), let us substitute equa- 
tions (W54), (W55) and (W14) into (W35). Now we see, that ipit) = 22)i2(t), (j){t) = ^22{t), i.e. 
we get a complete coincidence between our GFPE p.26p and Wang's GFPE (W29). 
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The PDF is evaluated directly through relations (|2TT8)) and ((27T21) with Kj = 0: 



fix,v,t) = 



: exp 



47r^OTii(<)OT22(<)-9?ti2(i)2 \4r!2 (Win WWlaaW - 9Ki2(i)2) 

X I 97lii (i) e''^/^vn - vo cosh (m/2) D, + (2x00;^ + i;o7) sinh {nt/2) 

r -1 2 

+ 2^22 (i) -e'^*/^xr2 + xo cosh {nt/2) il + (2wo + xqj) sinh (fii/2) 
+ 29Jti2(i) e'^^/'^vn - Vo cosh (m/2) O + (2a;oa;2 + uo7) sinh {fit/ 2) 

-e'^/'^xn + xo cosh {nt/2) ft + (2-i;o + xqj) sinh (r2i/2) 



(3.30) 



where dJlij are given with equations (| 3. 2 314X^5)1 . 

4. Transition to Einstein-Smoluchowsl(i equation 



Now let us prove that the system considered in section 13.31 may be described with Einstein- 
Smoluchowski equation at high viscosity levels and at long times. 

When uj/j <C 1, we can neglect the time derivative of velocity and, therefore, the pair of 
Langevin equations p.l9p transforms into a single overdamped Langevin equation: 



dx LiP' , . 1,,/ N 

dr 7 7 



(4.1) 



The GFPE for such a system, according to equations (14.11) and (I2.16P has the following form: 
with 



^^" 7^+1 Iv r(l-/3) )■ 

Executing the same calculations as in the previous section, for the PDF we unfold: 



p{x,t) = 



2^j-Km{t) 



exp 



X — XqC" 



-t"V7 



where 
9Jl(i) 



1 



ct^ 



2^ 27a;2r(l-/3) 



E, 



7 



4an(i) 



—M 1-/3,2-/3,— 

1-/3 V 7 



(4.3) 



(4.4) 



(4.5) 



Similar results were obtained by M. Caceres in [15] for a stationary case [see equations (2.14) 
with (2.17) of the mentioned paper]. 

Expanding the coefficient (14.51) into a series at large (7^) 's and substituting it to the PDF (|4.4p 
yield: 



p(x,t) 



n^P/2^2~p I ^-p^i3^-2p^-t^-/^ 



V27rc 
X exp 



r(i-/3) ) 



xxQ-i^uj^-^f^e-'^ h 



2c 



cT(l - /3) 



• (4.6) 
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Now we return to the PDF for the most general case (J3.30I) and also expand it into a series at 
w/7 < K 7t: 



f{x,v,t) « Aexp 



^-l3^2-2P^-u.hh i-^e-t7^2-2/3 ^-/3 (^2 _ 3^2^ " 



\ cr(l-/3) cr(l-/3) 2c / 



\ cr(i-;3)2 cr(i-/3) 

'uj^-f-^e-'^^^/'^ivo + 2:07) woe-*T72-'9 ' 



2c 



c c 

' ^0-2 A-213 ~2i^^th „-t7^/9-2, ,4-2,3 \ 



Wo 



^f. ^^, ^^g 



e -'7^' "cj 



2c 



Xo 



2 f 7^t^4-2/3g-2.^^t/7 e-t7^/3-2^6-2/3 

2^; c 



WqXo 



^P-l^i-2l3Q-2u,^t/-i g_t^^/3_i^4-2/3 



2nV:^/ ^~2/3^/3^-/3g-^^t/7 g-t7(^^)-2ff ' 

~ 47rc r^ r(i-;3) r(i-/?)2 



(4.7) 
(4.8) 



Then, integrating it by v in the range of (—00; 00) and neglecting the terms of the higher magnitude 
of smallness than exp {^—u'^t/^^ we get: 



f{x,t) ex exp 



2/7^^4-2,3 ^-^^2/3^4-4/3e-^^t/7 \ ^^g^/3^4-2/3g-t<.V7 



2c 



cr(i - 13) 



(4.9) 



which fully corresponds to the PDF (I4.6p . and, therefore, proves the fact that the considered system 
at large times and strong friction may be described with Einstein-Smoluchowski equation. 



5. GFPE for overdamped harmonic oscillator with constant drift 

As a final application example of the presented technique, let us study the PDF of the thermo- 
dynamical work w in the stochastic system which consists of a particle inside a harmonic potential 
moving with constant velocity v^, , U = {k/2)[x — X{t)]'^, X{t) — w+i, x (i) is the particle's coordi- 
nate. Our aim is to get the transient fluctuation relation for such a system, which will demonstrate 
large-deviation symmetry properties in the PDF, and compare it to the classical case. 

The work w is defined as follows: 



'{i) = / dX— = / di' — — = -kv, / dt'{x - vj). 



dX 



dt' dX 



(5.1) 



Introducing y (t) — x{t) — v^t, for the overdamped Langevin equation and the equation for the 
thermodynamical work w (t) we have: 



§ = --yit) + Yit)-v^ 
at T 

dw , , , 

— = -kv^y [t) , 



(5.2) 
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where r = m'^/k. Alternatively, if we consider the plane [y^w), the coefficient matrix a of the 
system will be 

" 1/r 
few* 



(5.3) 



Since y (t) = j/q exp(— t/r), j/o = ^o is the initial position of the particle, 

w (t) = wo + yokv^T f e"*/'^ - 1 j , 
Then the evolution matrix 



e-*/^ 

kv,T fe"*/^ - 1) 1 



For the diffusion coefficients 1)ij (t) we have: 



S)l2(t) 



ckv,.-fH-^ (t/T [(/3 - 1)E^ (t/T) - 1] + r(2 - (3) (t/rf) 



2r(2 - /?) 

©22 (t) = 0. 

The generalized Fokker-Planck equation in this case will have the form: 



df(y,w,t) nj \df , Of ^ d^f „^ ,, d'^f 
\T ' oy ow oy^ oyow 



dt 



(5.4) 

(5.5) 

(5.6) 
(5.7) 

(5.8) 



Now, considering an initial condition / (y, w, 0) — n5 {y — j/o) 5 (w — wq), when yo — 0^ wq = 
we get for the PDF: 



fiw,t) 



2V^V^Wi) 



exp ■ 



[kv^T^ {t/T + e-*/" - 1) + ^ 

4OT22(t) 



(5.9) 



where 



at2 



2r(3 - /?) 



2 - e-*/^M 



1,3-/3, 



-«/ 



^ (2 - e-*/^) 



M 



2-/3,3-/3, 



(5.10) 



which fully corresponds to the results obtained in [37|. 

After the relaxation stage has passed, at i ^ r we find for the transient fluctuation relation: 



f{w,t) 
f{-w,t) 



exp 



T{3-I3)wt^-^ 



cmj 



(5.11) 



Thus, the fluctuation relation for the system subjected to a coloured noise with the slowly decaying 
power-law correlation function differs from that for ordinary Brownian motion. As we stated above, 
the classical case limit is revealed at /3 ^ 1. 

6. Conclusions 

In this paper we suggested a consistent method for derivation of the generalized Fokker-Planck 
equation for linear multidimensional Gaussian non-Markovian systems. Taking the case of the 
Gaussian systems with slowly decaying power-law correlations, we obtained the following results: 
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• Firstly, we constructed the solution of generalized Fokkcr-Planck equation, the probability 
density function, without solving it directly. 

• We derived generalized Fokker-Planck equation for free motion and constructed the proba- 
bility density function for spatially homogeneous and inhomogeneous cases. 

• For the case of the motion in a harmonic potential, the generalized Fokker-Planck equation 
and the probability density function were also obtained, and the results were compared to 
those of the other authors. 

• We show the equivalence in description of generalized Brownian motion in a harmonic poten- 
tial with generalized Fokker-Planck equation and generalized Einstein- Smoluchowski equation 
at high viscosity levels and at long times. 

• Finally, we investigated the probability density function for thermodynamical work in the 
stochastic system which consists of a particle inside a uniformly moving harmonic potential 
underlining strong differences in transient fluctuation relations for the generalized Brownian 
motion and the ordinary Brownian motion cases. 
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Vaara/ibHeHe piBHJiHHfl OoKepa-fl/iaHKa ra Roro po3b'5I30k p,n9\ 

f\\H\\AHV\X HeMapKiBCbKI/IX TayCOBHX CMCTeM 

O.hO. C/iKDcapeHKO 

iHCTMTyx TeopeTMHHoT ctDiBMKM Im. 0.1. Axiesepa HHI_| X0TI, VKpaiHa, 61108 XapKis, ey/i. AKafleMiMHa, 1 

y LiiCi poSori MM nponoHyeMO nocjiiflOBHi/ifi niflxifl flo Bi/iBOfly ysarajibHeHoro pisHaHHa OoKepa-H/iaHKa (VP- 
On) p,nft laycoBwx HewapKiBCbKHX npoL;eciB is CTaL4ioHapHHMi/i npupomeHHaMH. U,e\A niflxifl p,03Bonife no- 
SyflyBaTM cJDyHKLiJKi po3nofli/iy (OP) npoL4ecy 5e3 noTpe6i/i 6e3nocepeflHbo posB'^ByBarn VPOH. Mm 3acTO- 
coByeMO L4eM Mexofl fljia 3Haxofl>KeHH!R VPOH ra 0P fl/in Bi/ibHoro y3arajibHeHoro 5poyHiBCbKoro pyxy ja 
y3arajibHeHoro SpoyHiscbKoro pyxy b noreHi^iajii flyin Bi/inaflKy cxeneHeBoT KopejiaLiitiHoT 4>yHKLiii' LuyMy. Mi/i 

flOBOflMMO, mo p03rjl51HyTi CMCTeMM MO>KyTb OnMCyBaTI/lCJI y paMKaX piBH!RHH!R EMHUJTeMHa-CMO/iyXOBCbKOrO 
3a yMOB CM/lbHoT B'513KOCTi Ta Be/IMKMX HaciB. TaKO>K Ml/1 HOpiBHKDeMO pe3yjlbTaTM i3 OTpMMaHMMM JHlUMMI/l 

aBTopaMM. HapeujTi, mm o5HMCJiKieMO 0P xepMOflMHaMinHoT po6oTM y cxoxacTMHHJti CMcreMi, mo CKJiafla- 
eTbca 3 HacTMHKM y rapMOHJHHOMy noTeHL4iajii, nKi/i£i pyxaexbca 3 nocxitiHOK) lUBHflKicTK), ra nepeBipaeMO 
4>/iyKTyaLii£iHy reopeMy fljia po6oTi/i y xaKil/i ci/icjeMi. 

K/iiOHOBi c/iOBa: p'lBHaHHa 0OKepa-n/iaHKa, laycosa ci/iCTSMa, neMapKiBCbKa ci/iCTema, TepMOfli/inaMi'-iHa 
poBoTa, nepexiflHe (pjiyKTyaLiiMHe cniBBiflHOLueHHSj 
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